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ABSTRACT 

Classical  dynamics  calculations  are  performed  for  0(^P)  +  F[20(’Ai)  collisions  from  2  to  10  km 
s'^  (4.1  to  101.3  kcal  mol  ’),  focusing  on  product  internal  energies.  Several  methods  are  used  to  produce 
ro-vihrationally  state -resolved  product  cross  sections  and  to  enforce  zero-point  maintenance  from  analysis 
of  the  classical  trajectories.  Two  potential  energy  surfaces  are  used:  1)  a  recently  developed  set  of  global 
reactive  surfaces  for  the  three  lowest  triplet  states  which  model  OH  formation,  H  elimination  to  make  H  -i- 
OOH,  0-atom  exchange,  and  collisional  excitation,  and  2)  a  non-reactive  surface  used  in  past  classical 
and  quantum  collision  studies.  Comparisons  to  these  previous  studies  suggest  that  for  H2O  vibrational 
excitation,  classical  dynamics  which  include  Gaussian  binning  procedures  and/or  selected  zero-point 
maintenance  algorithms  can  produce  results  which  approximate  quantum  scattering  cross  sections  fairly 
well.  Without  these  procedures,  the  classical  cross  sections  can  be  many  orders  of  magnitude  greater  than 
the  quantum  cross  sections  for  exciting  the  bending  vibration  of  H2O,  especially  near  threshold.  The 
classical  cross  section  over-estimate  is  due  to  energy  borrowing  from  stretching  modes  which  dip  below 
zero-point  values.  For  results  on  the  reactive  surfaces,  the  present  calculations  show  that  at  higher 
velocities  there  is  an  unusually  large  amount  of  product  internal  excitation.  For  OOH,  where  40%  of 
available  collision  energy  goes  into  internal  motion,  the  excited  product  vibrational  and  rotational  energy 
distributions  are  relatively  flat  and  values  of  the  OOH  rotational  angular  momentum  exceed  J=100.  Other 
product  channel  distributions  show  an  exponential  fall-off  with  energy  consistent  with  an  energy  gap  law. 
The  present  detailed  distributions  and  cross  sections  can  serve  as  a  guide  for  future  hyperthermal 
measurements  of  this  system. 
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1.  INTRODUCTION 


Collisions  of  0(^P)  atom  with  H20(^Ai)  have  recently  been  investigated  using  hyperthermal 
crossed  heams^  in  the  energy  range  ~50-  80  kcal  mol  ’  (~7-9  km  s  ’).  At  these  energies,  several  reactions 
occur: 

(1)  OH  formation:  0(^P)  +  H20(X,’Ai)  ^  OH(X,^n)  +  OH(X,^n) ,  AHo  =  15.8  kcal  mor’,^  AE^~18 
kcal  mol  ’  ^ 

(2)  H-elimination:  0(^P)  +  H20(X’Ai)  ^  OOH(X,^A”)  +  H(X^S)  AH,,  =  53.4  kcal  moP’  ^  AE^  ~60  kcal 
mof’  ^ 


(3)  0-atom  exchange:  0(^P)  -i-  H20’(X,’Ai)  ^  O’(^P)  +  H20(X,’Ai),  AH,,  =  0.0  kcal  mol  ’,  AE^  ~33  kcal 
mof’^ 

In  addition,  non -reactive  collisional  excitation  occurs: 

(4)  Collisional  excitation:  0(^P)  -i-  H20(X,’Ai,  (viV2V3,JKaKe))  ^  0(^P)  -i-  H20(X,’Ai(vi  V2  Va  J  Ka  Ke’)) 

In  Eqs.  (1-3),  the  AE^  values  are  approximate  threshold  energies  based  on  the  surfaces  of  Ref.  (3)  and 
include  zero-point  energy. 

Studies  of  0(^P)  -i-  H20(’Ai)  hyperthermal  collisions  go  back  nearly  40  years  to  measurements  of 
infrared  radiation  in  shock  tubes.®  In  the  1980’s,  cross  sections  for  vibrational  excitation  were  calculated 
using  quasi-classical  trajectory  (QCT)  methods^  on  a  global  non-reactive  potential  surface*.  These  cross 
sections  were  in  fair  agreement  with  the  limited  experimental  data  available.  Slightly  later,  approximate 
quantum  scattering  calculations  with  the  same  surface  were  performed.®  The  QCT  and  approximate 
quantum  cross  sections  were  orders  of  magnitude  apart,  even  at  relatively  large  collision  velocities.  These 
differences  have  made  it  difficult  to  interpret  the  measurements.  They  also  present  a  large  puzzle  in 
assessing  the  accuracy  of  classical  versus  quantum  methods. 

We  note  that  recently  there  have  been  several  advances  in  classical  dynamics  methods  for 
determining  state-resolved  excitation  and  reaction  cross  sections  in  polyatomic  systems.  These  have  led  to 
impressive  results  for  several  systems:  CH4  -1-  Cl  A  HCl  -1-  CHa,’”  OH  -1-  D2  ^  HOD  -1-  D,”  (H20)2  and 
(D20)2  dissociation,’^  and  rotational  excitation  of  H2O  by  H2.’*  In  particular  for  the  reaction,  OH  -1-  D2  ^ 
HOD  -I-  D,  Sierra  et  use  standard  QCT  methods,  plus  two  key  advances:  determination  of  vibrational 
quantum  numbers  through  analysis  of  product  momenta  and  coordinates,’"’  ’®  and  a  Gaussian  binning 
procedure  which  weights  trajectories  where  classically  determined  vibrational  quantum  numbers  are  close 

2 


Distribution  A:  Approved  for  pubiic  reiease;  distribution  uniimited. 


to  “allowed”  quantum  integer  values/^’’  With  these  procedures,  Sierra  et  obtained  excellent 
agreement  with  angle-velocity  measurements  of  product  D-atoms,  showing  distinct  vibrational  structure 
of  product  HOD  bending  and  stretching  modes.  For  (H20)2  and  (D20)2  dissociation,  Czako  et  al.  use 
similar  binning  procedures,  except  the  Gaussian  weights  are  determined  by  vibrational  energies  instead  of 
quantum  numbers.  QCT  methods,  which  reproduce  more  accurate  benchmark  quantum  scattering  results 
and  have  been  tested  against  experimental  measurements,  could  be  a  practical  alternative  to  quantum 
scattering  methods. 

The  primary  purpose  of  this  paper  is  to  present  results  of  0(^P)  -i-  H20(’Ai)  QCT  calculations 
using  these  recently  developed  classical  dynamics  methods.  The  focus  is  on  product  internal  energies  and 
ro-vibrational,  state-specific  cross  sections  for  Eqs.  (1-4).  Particular  attention  is  given  to  vibrational 
excitation  in  Eq.  4  to  better  understand  previous  hyperthermal  measurements  and  modeling.  In  one  set  of 
QCT  calculations,  the  work  of  Refs.  7  and  9  is  re-examined.  In  this  previous  work,  conventional 
histogram  binning  QCT^  and  approximate  quantum  scattering^  calculations  for  vibrational  excitation,  Eq. 
(4),  were  performed  using  the  non-reactive  0(^P)  -i-  H20(^Ai)  surface  of  Refs.  7,8,  and  18.  We  wish  to 
understand  the  more  than  order  of  magnitude  difference  in  vibrational  excitation  cross  sections  from  these 
two  studies  in  light  of  recent  advances  in  classical  dynamics  methods.  Eor  the  first  set  of  QCT 
calculations  reported  here,  the  same  non-reactive  surface  of  Ref.  7  was  used.  Ro-vibrational  H2O 
excitation  cross  sections,  Eq.  4,  were  determined  with  conventional  histogram  binning  and  Gaussian 
binning  at  collision  velocities  2,  4,  6,  and  8  km  s  ’  (4.1,  16.2,  36.5,  and  64.8  kcal  moP'),  with  10® 
trajectories  per  velocity.  Eor  these  QCT  calculations,  a  number  of  different  methods  were  used  and 
compared  for  analysis  of  coordinates  and  momenta  to  determine  product  vibrational  quantum  numbers.  A 
chief  finding  of  these  QCT  calculations  is  that  Gaussian  binning  dramatically  decreases  the  vibrational 
excitation  cross  section,  (Eq.  4),  for  the  H2O  bending  mode,  compared  to  conventional  histogram  binning, 
especially  near  threshold.  With  conventional  histogram  binning,  the  bending  mode  will  tend  to  borrow 
energy  from  stretching  vibrations  which  dip  below  their  zero  point  energy,  and  the  bending  excitation 
cross  sections  become  un-physically  large.  The  Gaussian  binned  QCT  results  were  in  fair  agreement  with 
the  approximate  quantum  scattering  calculations  of  Ref.  9  using  the  same  surface.  These  results  indicate 
that  the  several  orders  of  magnitude  difference  between  the  QCT  H2O  bending  cross  sections  of  Ref.  7 
and  the  approximate  quantum  scattering  cross  sections  of  Ref.  9  using  the  same  surface  are  due  to 
quantum  effects  that  can  be  effectively  accounted  for  with  Gaussian  binning.  Eor  exciting  the  stretching 
modes,  we  were  unable  to  reproduce  the  QCT  cross  sections  of  Ref.  7  with  the  same  QCT  histogram 
binning  procedures  and  other  parameters  used  in  that  study.  The  stretching  vibrational  excitation  cross 
sections  we  obtained  with  histogram  binning  were  much  smaller  in  magnitude  than  those  of  Ref.  7  and 
were  much  more  in  line  with  the  approximate  quantum  scattering  results  of  Ref.  9.  The  stretching 
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vibrational  excitation  cross  sections  changed  little  with  Gaussian  binning.  We  performed  approximately 
200  times  more  trajectories  per  velocity  than  Ref.  7,  but  otherwise  used  very  similar  methods  to  generate 
the  cross  sections. 

In  a  second  series  of  calculations,  the  recently  developed  0(^P)  +  H20(^Ai)  surfaces  of  Ref.  3 
were  used,  including  all  reactive  channels  (Eqs.  1-3),  for  the  three  lowest  triplet  electronic  states.  The 
surfaces  were  fits  of  CASSCF-I-MP2  electronic  structure  calculations  using  a  polynomial  function  of  the 
bond  distances,  with  built-in  permutational  symmetry  of  like-atoms.'®  These  calculations  were  performed 
at  velocities  of  2,  4,  6,  8,  and  10  km  s  *  (4.1,  16.2,  36.5,  64.8,  101.3  kcal  mol  '),  and  include  ~10^ 
trajectories  per  surface  per  velocity.  As  with  the  first  set  of  QCT  calculations,  Gaussian  binning 
dramatically  decreases  the  non-reactive  vibrational  excitation  cross  section,  (Eq.  4),  for  the  H2O  bending 
mode,  compared  to  conventional  histogram  binning,  especially  near  threshold.  In  contrast,  conventional 
histogram  binning  and  Gaussian  binning  gave  similar  results  for  the  reactive  channels,  Eqs.  (1-3),  except 
very  close  to  threshold.  Eor  OOEl,  where  40%  of  available  collision  energy  goes  into  internal  motion,  the 
excited  product  vibrational  and  rotational  energy  distributions  are  relatively  flat,  and  values  of  the 
rotational  angular  momentum  can  reach  values  greater  than  J=100.  Other  product  channel  distributions 
show  an  exponential  fall-off  with  energy  consistent  with  an  energy  gap  law.  We  also  compare  these 
results  to  the  first  set  of  QCT  calculations  using  the  surfaces  of  Refs.  7,  8,  and  18  to  examine  the 
sensitivity  of  the  cross  sections  to  the  potential  used.  For  excitation  of  the  bending  mode,  there  is  good 
agreement  between  the  two  QCT  Gaussian  binning  results,  except  very  close  to  threshold.  For  stretching 
excitation,  both  surfaces  give  small  cross  sections,  but  there  are  significant  differences.  This  is  mainly  due 
to  the  reactive  nature  of  the  Ref.  3  surfaces  which  break/form  bonds  as  the  0-H  bonds  change.  For  the 
two  sets  of  QCT  calculations  reported  here,  the  experimental  results  of  Ref.  6  fall  well  outside  the  quoted 
experimental  uncertainty  and  statistical  uncertainty  of  the  present  calculations.  As  a  further  test  of  the 
quality  of  the  potential  surfaces,  we  also  performed  QCT  rate  constant  calculations  at  fixed  temperatures 
for  OH  formation,  Eq.  (1)  and  the  reverse  of  Eq.  (1),  where  there  is  a  great  deal  of  measured  data.^°'^^  For 
these  rate  constants,  the  present  results  agreed  well  with  measurements. 

Since  H20('Ai)  zero-point  energy  (~12.9  kcal  moP')^"'  can  be  a  significant  fraction  of  the 
collision  energy,  especially  at  the  lower  collision  velocities  studied,  we  also  implemented  procedures  to 
maintain  zero-point  energy  during  the  classical  trajectory,  following  the  work  of  Refs.  25  and  26.  Zero- 
point  maintenance  in  QCT  calculations  may  be  important  in  certain  kinds  of  reactions,  especially 
involving  collision  complexes  with  long  lifetimes.  However,  the  implementation  of  zero-point 
maintenance  has  to  be  performed  with  great  care.  For  the  O  -1-  H2O  collisions  studied  here,  in  many 
instances,  we  found  that  it  is  difficult  to  cleanly  separate  vibrational  and  rotational  type  motion  during  the 
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collision.  When  the  O  and  H2O  reactants  were  close  together,  this  sometimes  resulted  in  ‘low  frequency 
vibrational  modes’  being  used  to  help  maintain  zero-point  vibrational  energy,  and  a  leaking  of  what  may 
be  rotational  energy  into  vibration.  As  a  result,  energetically  closed  vibrational  excitation  channels 
(channels  not  energetically  allowed  quantum  mechanically)  close  to  threshold  began  to  pick  up  significant 
probability.  When  we  restricted  the  zero-point  maintenance  procedures  to  portions  of  the  trajectory  where 
the  reactants  and  products  were  well-separated  and  there  was  better  separation  of  rotational  and 
vibrational  motion,  the  results  closely  matched  the  Gaussian  binning  results.  Our  conclusion  is  that  the 
zero-point  maintenance  procedures  when  applied  to  separate  product  and  reactant  portions  of  the 
trajectory  were  attempting  to  enforce  the  same  kind  of  partitioning  of  energy  among  the  vibrational 
modes  that  result  from  Gaussian  binning  procedures  for  the  processes  studied  here. 

The  paper  proceeds  as  follows.  Section  II  describes  details  of  the  QCT  calculations,  including  the 
methods  for  analysis  of  products  for  assignment  of  polyatomic  ro-vibrational  quantum  numbers.  Section 
II  also  describes  procedures  used  in  Gaussian  binning  and  zero-point  maintenance,  which  are  both  closely 
related  to  determination  of  the  polyatomic  quantum  numbers.  Section  III  presents  results  of  QCT 
calculations  on  the  non-reactive  surface  of  Refs.  8  and  18,  focusing  on  vibrational  excitation.  Results  are 
compared  to  the  QCT  calculations  of  Ref.  7  and  the  approximate  quantum  scattering  calculations  of  Ref. 
9.  Section  IV  presents  results  of  QCT  calculations  using  the  recently  developed  0(^P)  -1-  H20(*Ai) 
surfaces  of  Ref.  3,  focusing  on  the  internal  product  energies  of  processes  in  Eqs.  (1-4).  Section  V  gives 
conclusions. 

II.  METHODS 

A.  POTENTIAL  SUREACES 

Two  different  potential  energy  surfaces  were  used  for  the  0(^P)  -1-  H20('Ai)  dynamics 
calculations.  One  was  the  non-reactive  0(^P)  -1-  H20(’Ai)  surface  developed  by  Bartlett,  Redmon,  Schatz, 
and  Garrett  (BRSG).^’*’^*  This  surface  is  a  Sorbie -Murrell  many-body  fit  of  ab  initio  MBPT  calculations,^^ 
including  an  analytical  fit  to  the  non-bonding,  pair-wise  and  three  body  interactions  between  O  and  each 
atom  in  the  H2O  molecule  as  well  as  an  analytical  fit  to  the  ab  initio  surface  of  H2O.  The  authors  of  the 
BRSG  surface  note  that  it  is  best  suited  to  collision  energies  up  to  3.0  eV  (69.2  kcal  mol  ').’  As  discussed 
in  the  previous  section,  the  BRSG  surface  was  previously  used  for  quasi-classical  trajectory  calculations* 
and  quantum  scattering  calculations.^  While  there  have  been  some  modifications  to  the  BRSG  surface  by 
Schwenke  for  QCT  studies,’*  the  present  work  uses  the  original  BRSG  surface  to  enable  a  more  direct 
comparison  of  results  with  Refs.  7  and  9. 
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The  second  0(^P)  +  H20(’Ai)  potential  was  recently  developed  by  Conforti,  Braunstein,  Braams, 
and  Bowman  (CBBB).^  The  global  CBBB  surface  is  actually  a  set  of  the  three  lowest  triplet  electronic 
states  for  the  0(^P)  +  H20(’Ai)  system.  Approximately  10^  ab  initio  single  point  energies  were  calculated 
using  CASSCF+MP2  theory,  which  sampled  geometries  including  reactive  (Eqs.  1-3)  and  non-reactive 
(Eq.  4)  collisions.  In  addition,  many  configurations  were  included  for  high-energy,  short  bond  distance 
interactions  in  order  for  the  resulting  potential  energy  surface  to  be  applicable  at  high  collision  velocities. 
Eor  the  fit,  basis  functions  of  the  six  inter -nuclear  bond  distances  yielded  approximately  500  terms  for 
each  electronic  surface.  The  fitting  procedure  is  described  in  great  detail  in  Ref.  19.  The  CBBB  surfaces 
are  best  suited  for  collision  energies  up  to  approximately  150  kcal  moT^ 

Since  we  wish  to  understand  the  sensitivity  of  the  dynamics  simulations  to  the  potential  surfaces 
used,  it  is  appropriate  to  examine  differences  between  the  BRSG  and  CBBB  surfaces.  In  Eigure  1,  contour 
plots  are  shown  of  the  O  -i-  H2O  system  using  the  BRSG  (Eigure  la  and  Ic)  and  CBBB  surfaces  (Eigure 
lb  and  Id).  Eor  Eigure  la  and  lb,  the  E120  molecule  is  close  to  the  equilibrium  geometry,  and  the  0-0 
bond  distance  is  varied  within  H2O  plane.  Eor  Eigure  Ic  and  Id,  the  H2O  molecule  is  such  that  one  of  the 
0-H  bonds  is  stretched  from  its  equilibrium  value,  and  the  0-0  bond  distance  is  varied  in  the  H2O  plane. 
Eor  the  CBBB  surface  in  Eigure  Id,  the  lower  energy  region  indicates  the  formation  of  an  OH  bond.  The 
corresponding  BSRG  geometries  in  Eigure  Ic  show  a  deep  well  on  each  side  of  the  H-atom  that  has  been 
stretched.  The  BRSG  and  CBBB  surfaces  in  Eigure  la  and  Eigure  lb  very  closely  resemble  each  other. 
However,  the  surfaces  are  clearly  different  in  Eigure  Ic  and  Id.  We  note  that  the  BRSG  surface  is  not 
intended  to  be  reactive.  But  8  km  s'^  (64.8  kcal  moP')  collisions  should  access  these  geometries,  and  so 
the  resulting  dynamics  may  be  significantly  different. 

B.  Classical  dynamics  formulation 

In  this  section  we  review  several  procedures,  recently  developed  by  others  and  used  in  the  present 
work  for  computing  QCT  state -resolved  cross  sections  for  polyatomics.  We  start  with  methods  for 
separating  rotational  and  vibrational  motion  in  polyatomic  molecules.  This  separation  leads  to  algorithms 
for  determining  vibrational  quantum  numbers  and  to  approaches  for  maintaining  zero-point  energy  in 
classical  dynamics  calculations.  Einally,  we  review  Gaussian  binning  procedures  for  determining  state- 
resolved  cross  sections. 

1.  Separation  of  rotational  and  vibrational  motion 

Although  methods  for  determining  vibrational  quantum  numbers  for  diatomic  products  in  QCT 
calculations  are  well-defined,^^  corresponding  procedures  for  polyatomic  molecules  are  less  well-known. 
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This  is  mainly  due  to  the  fact  that  polyatomic  products  require  an  extra  step  of  separating  rotational  and 
vibrational  motion.  Several  authors have  discussed  procedures  for  this  separation.  In  the  present  work, 
we  follow  the  approach  of  Rhee  and  Kim,^°  a  general  formulation  which  leads  straightforwardly  to 
algorithms  for  zero-point  maintenance  and  to  several  alternate  determinations  of  vibrational  quantum 
numbers. 

The  Rhee  and  Kim  approach  (Ref.  30)  starts  by  defining  mass  scaled  Cartesian  coordinates  for  a 
system  of  N  atoms  in  the  center  of  mass  frame  with  the  center  of  mass  velocity  removed.  We  denote  the 
atomic  coordinates: 


X  3  [Xa^m^,  a  =  1,2, ...  ,N,i  =  x,y,  z],  (5) 

where  is  the  atomic  mass  of  atom  a.  For  convenience,  vector  notation  is  used.  Therefore  the  symbol, 
X,  denotes  a  vector  of  three  Cartesian,  space  fixed  coordinates  for  each  of  the  atoms,  a.  We  define  the 
angular  velocity,  Ui  ,  in  the  usual  way: 

"  =  Za  X  x'a  (6) 

where  the  inertia  tensor,  7,  is  defined  as: 

^  ~  Za[l  (J^a  ■  ^a)  ']  •  (7) 

The  separation  of  vibrational  and  rotational  motion  following  Rhee  and  Kim^°  starts  by  writing  the 
Cartesian  velocity  in  terms  of  pure  vibrational  and  rotational  components  that  can  be  separated: 

Xq.  -f  Uq.  -f  (l)j-  X  Xq..  (8) 

where  the  sum,  -f  Uq.,  is  the  pure  vibrational  component  and  ciy  X  x^  is  the  pure  rotational 
component.  The  vector,  u^,  is  defined: 

Ua  =  dir-)  X  Xa-  (9) 

The  vector,  dir-,  is  defined  in  terms  of  the  dynamic  equilibrium  vectors,  , 

dir  =  ■  Za  ««  X  Xa  (10) 

where, 

7  =  I.a[^  (.Xa  ■  da)  —  Xada  ■]  .  (11) 

The  dynamic  equilibrium  vectors  are  extremely  important  to  the  formulation.  As  we  shall  see,  the 
dynamic  equilibrium  vectors  are  determined  by  applying  a  rotation  to  the  molecule’s  static  equilibrium 
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vectors  so  that  they  align  with  the  present  position  in  the  trajectory  of  the  rotating  molecule.  Defining  the 
vibrational  displacement  vector,  |q., 

(12) 

the  so-called  Eckart  condition  (Ref.  31),  which  is  really  a  restatement  of  angular  momentum  conservation, 
is: 

Y.a  aaX  L  =  0.  (13) 

The  dynamic  equilibrium  vectors,  ,  are  found  by  the  Eckart  transform, 

aa  =  Pc  a,  (14) 

where  f  is  a  matrix  which  transforms  the  fixed  molecular  equilibrium  Cartesian  positions,  c^,  to  the 
dynamic  equilibrium  positions,  The  Eckart  transform^'  in  the  above  equation  essentially  rotates  the 
molecular  equilibrium  to  align  with  the  present  position  of  the  molecule,  so  that  any  displacements  are 
purely  vibrational.  Thus  the  vibrational  motion  is  cleanly  isolated.  Once  dii.  is  found,  the  vectors  v„  and 
u„  can  then  be  solved  for  in  terms  of  known  vector  quantities  using  the  above  equations. 

2.  Determination  of  quantum  numbers 

Now  that  the  vibrational  motion  of  the  polyatomic  has  been  isolated  in  the  vectors  v„  +Ucc,  the  next 
step  is  to  separate  the  vibrational  motion  into  specific  modes  and  to  determine  energies  and  semi -classical 
quantum  numbers  of  these  modes.  This  is  done  by  diagonalizing  the  mass -weighted  second  derivative  of 
the  potential  matrix  (Hessian)  at  each  time-step.  The  eigenvectors,  Li,  are  used  to  determine  the 
momentum  (Q)  and  position  (Q)  of  each  vibrational  mode,  indexed  by  /: 

Qi=li-{v  +  u),  (15) 

Qi=Lr  I  (16) 

There  are  several  approximate  ways  the  vibrational  potential  energy  of  a  particular  mode  can  be  defined. 
In  the  harmonic  approximation  (HA),^°  the  mode  potential  energy  is: 

(17) 

In  the  elimination  method  (EM),^°  the  mode  potential  energy  is: 

Vvib.i  «  E  (a  +  QiLi)  -  Via).  (18) 

Eor  both  methods  the  normal  mode  kinetic  energy  is: 
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(19) 


T'vib.l  —^Ql- 

The  total  vibrational  energy  of  a  particular  mode  is: 

^vib,l  =  T'vib.l  +  ^vib,l-  (20) 

The  quantum  number  of  a  particular  mode  is  then  defined  as: 

Kib.i  =  (21) 

^mode 

where  £'mode  is  the  harmonic  energy  of  that  mode,  and  int  stands  for  integer.  In  practice,  we  estimate  the 
vibrational  mode  energy  by  averaging  the  energy  of  the  collided  product  over  several  vibrational  periods, 
and  then  we  assign  quantum  numbers  for  a  particular  trajectory  based  on  this  average.  During  the 
averaging,  the  summed  mode  energies  would  sometimes  differ  greatly  from  the  total  vibrational  energy, 
mainly  due  to  geometries  where  the  approximate  mode  potential  energies  were  poor.  We  therefore  used  a 
cutoff  energy  of  10  kcal  mol  '.  If  the  summed  mode  energies  and  total  energies  differed  by  more  than  this 
cutoff,  then  the  analysis  point  was  not  included  in  the  average. 

Schatz  and  coworkers  use  a  different  method  to  compute  quantum  numbers  from  the  classical 
eigenmode  momentum  and  position.'^  They  determine  the  classical  action  integral  involving  the 
momentum  and  position  of  a  particular  mode  over  several  vibrational  periods  and  then  relate  this  to  a 
quantum  number  following  the  WKB  method.  This  is  done  with  an  FFT  method.  In  our  work,  we  also 
compute  quantum  numbers  with  this  procedure,  denoted  as  “ACT”  cross  sections.  We  have  found  all 
methods  (HA,  EM,  and  ACT)  in  general  give  similar  results. 

3.  Zero-point  maintenance 

In  our  work,  we  have  also  employed  algorithms  to  maintain  zero-point  energy  during  and  after 
collisions  following  the  work  of  several  groups. These  algorithms  also  start  with  the  separation  of 
rotational  and  vibrational  motion.  We  start  by  defining  the  momentum  vector  for  the  system  which  is  the 
time  derivative  of  the  mass  scaled  Cartesian  coordinates: 

p  =  X.  (22) 

We  then  generate  the  projected  Hessian  matrix  as  discussed  in  Bonhommeau  and  Truhlar  (Eqs.  7  and  8).^® 
This  is  done  by  taking  the  second  derivative  potential  matrix  (Hessian),  creating  a  projection  operator, 
and  then  using  the  operator  to  project  out  the  translation  and  rotation  from  the  Hessian  to  isolate  the 
vibrational  motion.  We  then  calculate  the  eigen-frequencies  and  eigenvectors  of  the  projected  Hessian. 
The  eigen-frequencies  are  labeled,  H;,  and  the  eigenvectors  are  labeled. 
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,  for  the  l-th  vibrational  component.  We  proceed  by  defining  the  normal  mode  momenta,  P,  for  all 
the  real,  non-zero  vibrational  modes,  /. 


p,  =  iFr’ 


•p- 


(23) 


In  practice,  we  impose  a  lower  energy  limit  of  50  cm'^  to  a  vibrational  mode.  Therefore,  we  only  treat  the, 
3N-q,  (q>6),  modes  above  this  lower  limit.  The  instantaneous  energy  of  mode  I  is  defined  as: 


E,  =  0.5 


+  (: 


aiJ 


where  the  D;are  projected  mass-scaled,  first  derivatives  of  the  potential: 

Di  =  ■  W, 


(24) 


(25) 


where  in  vector  notation: 


W  3  { 


a 


1.2  =  x,y,z}. 


We  now  define  the  local  vibrational  mode  energy,  in  atomic  units,  as: 


(26) 


^zpe,l  —  0-6  fl;  ■  (27) 

We  use  the  following  “mTRAPZ”  criterion^^  for  when  to  enforce  zero-point  maintenance: 

(28) 

That  is,  whenever  the  sum  of  instantaneous  mode  energies  is  less  than  or  equal  to  the  sum  of  local 
vibrational  mode  energies,  we  re-scale  the  momenta.  We  do  this  in  the  following  way.  For  those  modes 
that  dip  below  their  zero-point  energy,  the  normal  mode  momenta  are  re-defined  to  compel  zero-point 
maintenance: 


p;  =  sgn(P,)>,  -  (^) 


(29) 


Total  energy  conservation  is  ensured  by  scaling  the  momenta  of  the  other  (non-violating)  modes: 

PI  =  PPu 

where. 


(30) 


l^  = 


^3N-q 

'-‘1=1 


o2  ■sT'nvio  r>f2 
I  Li=i  t'l 


3N-q 


■nvio+l^l 


(31) 
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and  where  nvio  is  the  number  of  zero-point  energy  violating  modes.  We  note  that  there  must  he  enough 
momenta  in  the  non-violating  modes  to  achieve  a  rescaling  (/?  must  he  real  and  greater  than  zero). 

The  last  steps  in  the  procedure  involve  transforming  hack  from  the  projected  space  {P’)  to  the 
Cartesian  space  {p  ’).  This  is  done  using  Eq.  29  of  McCormack  and  Lim,^^ 

Pa  =  ("a  X  ^a)  +  Pi  («),  (32) 

where  uIq.  is  the  angular  velocity  of  atom  a.  Finally,  we  note  that  due  to  numerical  round-off  and  putting  a 
non-zero  energy  limit  on  the  allowed  vibrational  modes,  there  will  be  some  small  energy  error  due  to  this 
final  transformation.  So  we  make  a  final  scaling  of  the  Cartesian  momenta,  defining  the  scaling  factor,  y. 


Y  = 


(33) 


where  KE  is  the  total  kinetic  energy  before  and  KE  ’  is  the  total  kinetic  energy  after  enforcing  zero-point 
maintenance.  So  that  the  final  Cartesian  momenta  are: 


Pa=YPa-  (34) 

This  final  scaling  is  usually  very  small,  perhaps  one  part  in  10^.  However,  with  trajectories  with  tens  of 
thousands  of  time  steps,  these  small  errors  can  mount  up  and  become  large.  To  further  ensure  that  the 
zero-point  algorithm  does  not  perturb  the  system  too  much,  we  restrict  the  value  of  yto  be  between  0.99 
and  1.01  in  a  single  time  step. 

4.  Histogram  and  Gaussian  binning 


We  follow  the  approach  discussed  in  Czako  and  Bowman^^  and  Bonnet  and  Rayez,'®  and  extend  it 
slightly  to  incorporate  reactive  collisions  where  several  product  channels  are  possible.  The  Gaussian 
binning  approach  starts  with  the  expression  for  the  vibrationally  resolved  cross  section  to  state  v: 

(Ty  T^t^maxPrPv)  (35) 


where  b^ax  is  the  maximum  impact  parameter.  In  the  above  equation. 


Pr  = 


N 


traj 


(36) 


is  the  probability  of  particular  reaction  channel,  where  W  is  the  number  of  trajectories  which  correspond 
to  a  particular  reaction  channel  (such  as  collisional  excitation)  and  Ntraj  is  the  total  number  of  trajectories. 
For  conventional  ‘histogram  binning’  (HB),  P„is  defined, 

Pv™  =  %  (37) 
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where  Ny  is  the  number  of  trajectories  whose  product  ends  up  in  particular  vibrational  state  v.  This  gives 
the  regular  expression  for  the  conventional  histogram  binning  (HB)  cross  section: 


HB 


=  rthl,, 


Ny 


N. 


traj 


(38) 


where  the  ‘weight’  of  each  trajectory  is  equal  to  1.0.  The  Gaussian  binning  approach  redefines  the 
probability  of  a  particular  state  in  terms  of  Gaussian  weights.  So  that 


pGB  _ 


ZpfiGp(v) 

2fiiZp=iGp(v)’ 


(39) 


where 


G(v)  = 


2  (E(v0-^  2 
^  2E(0)  ^ 


(40) 


In  the  above  equations,  Ny  is  the  number  of  trajectories  going  to  a  vibrational  product  state  v,  Ns  is  the 
number  of  product  vibrational  states,  E(  v’)  corresponds  to  the  energy  of  the  vibrational  state  at  the  end  of 
a  trajectory  and  E(  v)  corresponds  to  that  energy  where  the  quantum  number  is  exactly  an  integer  value. 
So  the  quantity  E( v)-E( v)  corresponds  to  the  ‘distance’  away  from  the  true  integer  quantum  state.  The 
parameter  /?  is  the  width  of  the  Gaussian  weight,  and  is  set  at  a  value  of  16.651  so  that  the  width  of  the 
weighting  function  is  0.1  out  of  1.0  (FWHM). 

5.  Rotational  quantum  numbers  and  diatomic  vibration 

For  the  rotational  quantum  numbers  of  diatomic  and  polyatomic  products,  J,  we  use  the  convention: 

;  =  l/p|-i/2,  (41) 

where  Jy,  is  the  angular  momentum  vector.  In  principle,  it  may  be  possible  to  determine  (JKaKc)  resolved 
cross  sections  for  1420  products,  as  shown  in  the  work  of  Ref.  13.  In  the  present  work,  we  do  not  follow 
the  products  long  enough  in  time  to  make  a  meaningful  estimate  of  the  Ka  and  Kc  quantum  numbers.  For 
determining  OH  diatomic  vibrational  quantum  numbers  we  evaluate  the  integral  for  the  semi-classical 
action,  following  standard  procedures. 

C.  Classical  dynamics  calculations 


QCT  calculations  using  a  classical  dynamics  code  written  by  the  authors  were  performed  using 
the  BRSG  surface  and  the  three  lowest  triplet  CBBB  surfaces.  In  the  calculations,  the  initial  relative 
orientation  of  the  O  -i-  H2O  complex  was  randomly  chosen  according  to  Ref.  33  with  an  initial  separation 
between  H2O  and  O  of  17.0  au.  and  a  maximum  impact  parameter  of  7.56  au.  For  these  calculations,  the 
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H2O  molecule  is  in  the  ground  vibrational  state  with  only  the  zero  point  normal  mode  vibrational  energy 
and  no  rotational  energy.  The  simulation  time  step  was  IxlO  ’®  s  and  a  velocity-verlet  integrator  was  used. 
Collision  energies  included  2,  4,  6,  and  8  km  s  ’  for  the  BRSG  surface  and  2,  4,  6,  8,  and  10  km  s  ’  for  the 
CBBB  surfaces.  Products  were  declared  when  bond  distances  and  angles  matched  equilibrium  values 
within  a  specified  criterion  and  were  separated  from  each  other  by  more  than  17.0  au.  Trajectories  were 
then  continued  and  product  analysis  occurred  for  5.12x10'’'’  additional  seconds.  Using  these  procedures, 
products  of  all  trajectories  were  successfully  assigned  for  the  CBBB  surfaces  at  8  km  s  ’.  For  the  BRSG 
surface  at  8  km  s  ’,  about  0.1  %  of  the  trajectories  failed  by  either  very  poor  energy  conservation  and/or 
dissociating  to  products  other  than  O  +  H2O,  while  at  2  km  s  ’  all  trajectories  on  the  BRSG  surface  were 
successful.  It  is  known  that  the  BRSG  surface  starts  to  fail  at  energies  greater  than  8  km  s  ’  which  could 
likely  be  due  to  the  surface  deficiency  shown  in  Figure  Ic,  and  we  may  be  accessing  some  of  those 
geometries.  Average  energy  conservation  defined  as  the  average  of  the  energy  change  magnitude  between 
the  start  and  stop  of  a  trajectory  at  8  km  s  ’  was  0.19  kcal  mol  ’  for  the  BRSG  surface  and  6.8x10'^  kcal 
mol  ’  for  the  lowest  CBBB  surface.  For  the  CBBB  surfaces,  cross  section  results  at  a  single  velocity  for 
each  electronic  state  were  summed  and  divided  by  three.  Table  I  lists  these  and  other  simulation 
parameters.  Table  II  lists  vibrational  frequencies  and  equilibrium  geometries  of  reactants  and  products  for 
each  surface  and  literature  values. 

Simulations  were  run  both  with  and  without  zero  point  maintenance,  as  described  earlier.  The 
zero-point  maintenance  simulations  were  run  in  two  ways:  with  zero-point  maintenance  for  the  product 
analysis  only  (ZPl),  and  with  zero-point  maintenance  enforced  during  the  entire  trajectory  (ZP2).  We 
defined  an  energy  lower  limit  of  50  cm  ’;  if  a  vibrational  mode  had  less  than  this  energy,  it  was  not 
assigned  as  a  vibration.  For  ZP2  simulations,  it  was  necessary  to  define  a  force  threshold.  If  the  Hessian 
matrix  elements  between  atoms  exceeded  this  threshold,  they  were  considered  part  of  the  same  molecule 
when  zero-point  procedures  were  applied.  Product  analysis  and  zero-point  maintenance  required  the 
calculation  of  a  matrix  of  potential  energy  second  derivatives  (Hessian).  Depending  on  the  complexity  of 
the  potential,  calculation  of  the  Hessian  could  drive  needed  computer  resources.  At  8  km  s  ’,  it  was  found 
that  without  zero  point  maintenance,  calculation  of  the  Hessian  for  product  analysis  alone  could  take  as 
long  as  the  trajectory  itself.  Zero-point  maintenance  through  the  whole  trajectory  (ZP2)  could  increase  the 
needed  computer  resources  by  two  times  or  more  compared  to  ZPl  calculations.  It  was  found  that 
enforcing  zero-point  maintenance  for  separated  products  only  versus  separated  reactants  and  products  did 
not  change  results,  and  saved  considerable  computer  time. 
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III.  QCT  Results  Using  the  BRSG  Surface 


In  this  section,  results  of  QCT  calculations  using  the  BRSG  surface  are  discussed,  focusing  on 
collisional  excitation  of  H2O  fundamental  vihrational  modes  (Eq.  4).  Figure  2  shows  the  present 
calculated  cross  sections  versus  velocity  along  with  previous  calculations^’^  using  this  surface  and  a 
measurement  at  4.3  km  s"^  for  the  bending  and  asymmetric  stretch.®  The  present  results,  computed  at  2,  4, 
6,  and  8  km  s‘^  collision  velocity  (4.1,  16.2,  36.5,  and  64.8  kcal  mol  '  collision  energy),  are  shown  as  three 
lines  with  symbols  within  each  panel.  Each  colored  line  uses  a  different  method  to  find  the  vibrational 
quantum  number:  “QCT  ACT”  (black)  uses  the  Fourier  transform  method  of  Schatz  and  co-workers  to 
find  the  vibrational  action,'®  “QCT  HA”  (red)  uses  the  harmonic  approximation  of  Rhee  and  Kim,®"  and 
“QCT  EM”  (blue)  uses  the  elimination  method  of  Rhee  and  Kim.®"  The  top  row  of  four  panels  is  for 
exciting  the  fundamental  bend.  The  middle  row  panels  are  for  exciting  the  fundamental  symmetric  stretch. 
The  bottom  row  panels  are  for  exciting  the  fundamental  asymmetric  stretch.  For  a  given  row,  each  panel 
shows  the  present  results  for  a  different  method  of  binning  or  zero-point  maintenance:  “ZPO”  is  the 
baseline  QCT  method  and  refers  to  conventional  histogram  binning  and  no  zero -point  maintenance 
algorithm  employed,  “ZPl”  uses  conventional  histogram  binning  and  enforces  zero-point  maintenance  in 
the  product  analysis  step,  “ZP2”  uses  conventional  histogram  binning  and  enforces  zero-point 
maintenance  during  the  entire  trajectory,  and  “GB”  uses  Gaussian  binning  and  no  zero-point  maintenance. 
The  statistical  error  for  these  calculations  is  less  than  1%  for  cross  section  values  above  5.xl0'®'  m®  and 
about  10%  for  cross  sections  at  5.xl0'  m  .  However,  by  5.xl0‘  m  the  statistical  error  is  about  30  %. 
Results  at  or  below  this  cross  section  value  are  therefore  not  meaningful.  Furthermore,  the  Gaussian 
binned  results  suppress  many  trajectories,  so  the  effective  uncertainty  is  greater  than  described  above.  The 
two  results  labeled  “QCT”  (purple)  are  from  Ref.  7,  which  use  the  same  approximations  as  the  baseline 
histogram  binning  and  no  zero-point  algorithm  case  (“ZPO”  “QCT  ACT”)  and  the  same  surface,  but  are 
run  with  5,000  trajectories,  an  initial  rotational  energy  of  0.08  eV  =  1.8  kcal  mol  '  ,  and  a  maximum 
impact  parameter  of  5.3  au.  For  the  QCT  results  of  Ref.  7,  the  purple  line  with  symbols  represents 
tabulated  results,  and  the  purple  line  with  no  symbols  shows  results  obtained  by  fitting  to  a  functional 
form  as  shown  in  Figures  1  and  2  of  Ref.  7.  The  quoted  statistical  error  is  about  5-10  %  for  the 
fundamental  bend  mode  and  about  20-25  %  for  the  fundamental  stretch  modes.  The  results  labeled  VCC- 
lOS  (green  line  with  symbols)  are  from  approximate  quantum  scattering  calculations  of  Johnson  using  the 
same  BSRG  surface." 

We  first  discuss  the  present  QCT  results  for  the  bend  excitation  (top  row  of  Figure  2)  and 
evaluate  the  zero-point  maintenance  algorithms  and  binning  methods.  The  chief  result  is  that  there  is  a 
dramatic  lowering  of  the  cross  section  especially  at  lower  velocities,  over  two  orders  of  magnitude  at  4 
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km  s  ’,  in  going  from  the  ZPO  to  the  ZPl results  and  from  the  ZPO  to  the  GB  results.  The  ZPl  and  GB 
results  generally  match  each  other  well.  Examination  of  individual  trajectories  leading  to  bending 
excitation  shows  that  the  large  bending  ZPO  cross  sections  at  low  velocities  mainly  result  from  the 
borrowing  of  stretch  mode  vibrational  energy  which  goes  below  zero-point  values.  Maintaining  zero- 
point  energy  in  the  product  region  with  ZPl  mostly  corrects  this  problem.  The  GB  results  also  suppress 
trajectories  which  borrow  energy  from  stretching  modes  that  go  below  zero-point.  The  ZP2  cross  section 
matches  the  ZPl  and  GB  results  fairly  well  at  6  km  s  ’  and  8  km  s  ’,  but  remains  un-physically  large  near 
threshold.  (We  note  that  at  2  km  s  ’,  there  is  not  enough  collision  energy  to  allow  excitation  of  the 
bending  mode,  without  borrowing  from  zero-point  stretches.  Non-zero  cross  sections  for  bending 
excitation  at  these  velocities  are  good  indicators  of  a  deficiency  with  the  QCT  method  used.)  For  ZP2 
trajectories  leading  to  bending  excitation  at  lower  velocities,  it  was  found  that  in  the  strong  interaction 
region  where  the  O  and  H2O  are  near  each  other,  very  low  frequency  vibrations  of  the  O  -1-  H2O  complex 
were  highly  excited.  Before  the  products  entered  the  O  -1-  H2O  product  region,  the  low  frequency  modes 
would  transfer  energy  to  the  bending  vibration  leading  to  bending  excitation  for  a  few  trajectories,  enough 
to  boost  the  cross  section  at  2  km  s  ’  and  4  km  s  ’  to  relatively  large  values.  The  interpretation  of  this  was 
that  internal  rotation  motion  of  the  O  -1-  H2O  complex  was  being  mistakenly  changed  into  bending 
vibration  of  the  H2O  through  our  particular  implementation  of  the  zero-point  maintenance  algorithm  with 
ZP2  procedures.  It  may  be  possible  to  exclude  these  low  frequency  modes  by  setting  a  higher  energy 
cutoff  for  counting  vibrational  modes.  However,  it  was  difficult  in  practice  to  implement  a  satisfactory 
cutoff  for  the  entire  trajectory  which  did  not  exclude  what  was  obviously  the  H2O  bend  vibration. 

There  are  systematic  differences  between  the  QCT  HA,  QCT  EM,  and  QCT  ACT  results. 
Although  the  QCT  EM  mode  energies  at  any  given  time  step  are  more  accurate  than  the  QCT  HA  mode 
energies,  the  QCT  EM  mode  energies  tend  to  undershoot  the  true  mode  energies  for  the  fundamental 
modes.  This  behavior  follows  from  the  forms  of  the  EM  and  HA  potential  energies.  The  errors  for  the 
QCT  HA  mode  energy  were  generally  unbiased,  so  the  QCT  HA  averaged  estimates  are  probably  more 
accurate  than  the  QCT  EM  results,  at  least  for  the  fundamental  modes.  In  general,  the  QCT  HA  and  QCT 
ACT  results  agreed  well  for  cross  section  values  with  a  small  statistical  uncertainty.  For  these  reasons  we 
will  therefore  focus  mainly  on  the  QCT  ACT  results  with  ZPl  or  GB,  noting  cases  where  they  differ  from 
each  other. 

In  principle,  the  present  QCT  ACT  results  with  ZPO  (conventional  histogram  binning)  use  the 
same  methods  as  Redmon  et  alJ  and  the  results  should  be  directly  comparable.  The  present  results  for 
ZPO  (conventional  histogram  binning)  only  match  the  results  of  Redmon  et  al.^  at  8  km  s  ’.  The  present 
ZPO  results  are  smaller  than  those  of  Redmon  et  alJ  at  lower  velocities  by  up  to  a  factor  of  4  or  5  for  the 
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bending  mode,  much  larger  than  statistical  uncertainty  would  allow.  Indeed,  for  excitation  of  the 
stretching  modes  these  differences  are  much  larger.  We  attempted  to  exactly  duplicate  the  Ref.  7 
calculations  by  running  5,000  trajectories  at  4  km  s'^  with  the  same  impact  parameter  and  initial  rotational 
energy  following  Ref.  7.  The  resulting  fundamental  bend  cross  section  went  from  4.1x10'^’  m^  using  our 
original  parameters  in  Table  I  to  2.4x10'^^  m^  with  the  parameters  (maximum  impact  parameter  and  initial 
rotational  energy)  of  Ref.  7.  Although  the  new  value  agrees  even  less  well  with  Ref.  7,  this  difference  is 
on  the  order  of  the  statistical  uncertainty  of  the  5,000  trajectory  run.  It  is  possible  that  our  reconstruction 
of  the  BSRG  surface  is  in  error.  Other  than  this  possibility,  we  are  at  a  loss  to  explain  the  sometimes  large 
differences  between  the  present  ZPO  ACT  results  and  those  of  Ref.  7,  especially  for  the  stretching 
excitation  cross  sections. 

On  other  hand,  the  QCT  ACT  ZPl  and  QCT  ACT  GB  cross  sections  agree  well  with  the 
approximate  quantum  scattering  cross  sections  of  Ref.  9  using  the  same  BSRG  surface.  This  level  of 
agreement  between  QCT  GB  and  quantum  scattering  cross  sections  has  been  seen  in  other  systems. For 
O  +  H2O  bending  excitation,  the  GB  and  ZPl  approaches  properly  correct  for  an  extremely  large 
deficiency  in  conventional  QCT  methods.  We  have  attempted  to  optimize  the  value  of  the  weighting 
width  parameter,  /?,  used  in  the  GB  calculations  to  better  agree  with  the  quantum  scattering  results  of  Ref. 
9.  Better  agreement  can  be  achieved  at  velocities  of  4  km  s  ’  and  higher.  However,  this  improvement 
comes  at  the  expense  of  increased  cross  sections  at  2  km  s  ’,  which  are  not  strictly  allowed.  Finally,  we 
note  that  the  results  of  Dunn  et  al.^  are  more  than  10  times  greater  than  the  quantum  scattering  results  of 
Ref.  9  and  the  present  QCT  results  GB  and  ZPl  results. 

Figure  2  also  shows  cross  sections  for  exciting  the  fundamental  symmetric  stretch  and 
asymmetric  stretch  modes.  The  present  QCT  stretching  cross  sections  at  all  levels  of  theory  are  much 
smaller  than  the  present  bend  cross  sections.  For  the  stretches,  the  difference  between  conventional  QCT 
ZPO  and  QCT  ZPl/GB  results  are  much  smaller  than  for  exciting  the  bend.  For  the  asymmetric  stretch, 
the  QCT  EM  result  lies  about  an  order  of  magnitude  below  the  QCT  HA  and  QCT  GB  results,  and  we 
believe  this  is  due  to  the  systematic  underestimate  of  the  QCT  EM  mode  energy,  mentioned  earlier.  The 
present  QCT  ACT  and  QCT  HA  results  with  ZPl  and  GBO  are  within  a  factor  of  two  or  three  of  the 
approximate  quantum  scattering  results  using  the  same  surface  of  Ref.  9  near  8  km  s  ’.  At  6  km  s  ’  the 
present  results  are  zero  within  the  statistical  error,  and  clearly  more  trajectories  are  needed  to  make  a 
meaningful  comparison  at  this  velocity.  The  difference  between  the  QCT  results  of  Ref.  7  using  the  same 
surface  and  the  present  results  is  large,  amounting  to  more  than  an  order  of  magnitude  at  8  km  s  ’,  and  by 
inference  several  orders  of  magnitude  at  lower  velocities.  However,  at  4.3  km  s  ’,  the  results  of  Dunn  et 
al.^  for  the  asymmetric  stretch  have  fair  agreement  with  the  QCT  results  of  Ref.  7.  Taken  together,  our 
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impression  is  that  the  consistently  good  agreement  between  the  present  QCT  ZPl/GB  results  and  the 
approximate  quantum  scattering  results  for  all  fundamental  modes  indicates  that  our  implementation  of 
the  BSRG  surface  is  likely  correct,  and  that  either  a  QCT  ZPl  or  GB  approach  can  correct  for  deficits  in 
conventional  QCT  methods. 

To  help  understand  the  present  QCT  approaches  at  a  deeper  level,  Figure  3  shows  contours  of  the 
prohahility  for  exciting  hend  and  stretch  modes  as  a  function  of  the  classical  vihrational  action  for  the 
ZPO,  ZPl,  and  ZP2  approaches  at  2,  4,  and  8  km  s  ’  collision  velocities.  The  loglO  of  the  prohahility 
contour  values,  normalized  so  that  each  suh-panel  sums  to  1.0,  are  shown  over  a  very  large  dynamic  range 
to  highlight  low  value  cross  sections.  The  symmetric  and  asymmetric  stretching  actions  have  been 
averaged  to  a  single  effective  stretch  mode  action  for  ease  of  viewing.  The  ZPO  panel  results  show  that 
the  probability  of  exciting  the  bending  mode  (actions  >  0.5)  derives  mainly  from  trajectories  which  have 
borrowed  energy  from  the  stretch  modes,  which  dip  below  actions  equal  0.0.  From  these  plots,  it  is  clear 
how  the  ZPl  algorithm  reduces  the  amount  of  bending  excitation.  Gaussian  binning,  which  effectively 
multiplies  these  probabilities  by  peaked  functions  at  integer  values  of  the  actions,  achieves  a  similar 
effect.  It  is  also  clear  that  at  4  km  s  '  and  to  some  extent  at  2  km  s  ’,  the  ZP2  results  start  to  allow  the 
borrowing  of  stretch  mode  zero-point  to  excite  the  bend,  which  had  been  suppressed  in  the  ZPl  results. 

Figure  4  gives  a  different  view  of  the  vibrational  cross  sections.  Figure  4a  gives  vibrational 
energy  distributions  binned  at  2  kcal  moP’  for  the  4  collision  velocities  studied  for  ZPl.  These 
distributions  are  peaked  near  the  zero-point  energy  of  H2O  at  13.6  kcal  mol  ’  (for  the  BRSG  potential)  and 
fall  off  rapidly  with  energy.  Even  with  zero-point  maintenance  (and  Gaussian  binning),  there  are  many 
trajectories  which  fall  below  the  H2O  zero  point  energy,  giving  rise  to  non-zero  cross  sections  below  the 
peak  at  13.6  kcal  mol  ’.  As  the  collision  velocity  increases,  the  width  of  the  distribution  around  the  zero 
point  increases,  especially  at  8  km  s  ’  where  there  is  a  high  energy  tail  extending  up  to  30  kcal  mol  ’.  The 
increasing  width  of  the  distributions  indicates  increasing  vibrational  excitation  with  velocity,  but  gives 
little  quantitative  information.  Figure  4b  shows  vibrationally  state -resolved  cross  sections  at  6  and  8  km 
s  ’  versus  vibrational  energy.  At  8  km  s  ’,  there  appears  to  be  an  exponential  fall-off  of  the  cross  section 
with  vibrational  energy,  which  is  consistent  with  an  energy  gap  law.  Figure  4c  shows  the  angular 
distribution  cross  sections  for  the  ground  state  and  fundamental  vibrationally  excited  modes.  For 
vibrational  excitation,  the  H2O  is  clearly  back-scattered,  which  is  consistent  with  a  picture  of  a  hard 
collision. 

Figure  5  gives  an  overview  of  the  calculated  rotational  distributions  for  O  -1-  H2O  collisions.  Figure  5a 
gives  the  cross  section  of  low-lying  H2O  J  levels  summed  over  all  vibrational  levels  over  a  wide  energy 
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range.  As  indicated,  below  ~4  kcal  mol'^  results  from  Schwenke^*  are  shown  for  comparison.  Schwenke 
used  a  slightly  modified  BSRG  potential  and  did  not  include  H2O  vibration.  There  are  three  results  shown 
for  the  Schwenke  study:  standard  classical  dynamics  (symbols  without  lines),  full  quantum  scattering 
(symbols  with  solid  lines),  and  approximate  quantum  scattering  (symbols  with  dotted  lines).  The  three 
Schwenke  approaches  are  very  close  and  cannot  be  easily  distinguished  on  the  scale  of  the  figure,  except 
for  the  J=1  level,  where  the  standard  classical  results  overestimate  the  quantum  cross  sections  by  a  little 
over  a  factor  of  two.  Although  they  do  not  overlap  in  energy,  the  present  results  are  consistent  with  the 
corresponding  QCT  calculations  of  Schwenke.  Figure  5b  shows  calculated  J-resolved  cross  sections  for 
QCT  ACT  and  Gaussian  binning  at  the  4  collision  velocities  studied,  summed  over  all  vibrational  levels. 
At  the  highest  collision  energy,  J-levels  of  30  and  above  are  populated.  Figure  5c  shows  ro-vibrationally 
resolved  cross  sections  for  the  ground  and  fundamental  vibrational  levels  at  8  km  s  ’.  The  (000)  J-resolved 
cross  section  falls  off  fairly  rapidly  with  energy;  the  (010)  cross  section  is  much  flatter  versus  J-level, 
with  a  maximum  near  J=5  and  gentle  fall-off  until  there  is  a  cutoff  around  J=30.  The  statistics  for  the 
stretching  resolved  J-cross  sections  are  poor,  but  give  some  indication  that  they  are  relatively  flat  with  J- 
level.  The  vibrationally  resolved  rotational  energy  distributions  at  8  km  s  ’  are  shown  in  Figure  5d.  For 
any  single  rotational  energy  bin,  there  are  several  J  levels  which  contribute.  This  creates  a  rolling  average 
of  the  J-resolved  cross  sections  seen  in  Figure  5c.  For  the  (000)  and  (010)  levels,  except  for  the  lowest 
energy  bin,  the  distributions  fall  off  with  an  energy  gap  type  of  behavior.  For  the  stretching  modes,  again 
the  statistics  are  poor,  but  the  rotational  energy  distribution  is  flat. 

IV.  QCT  Results  Using  the  CBBB  surfaces 
A.  Ro-vibrational  excitation  of  H2O 

In  Figure  6,  cross  sections  for  excitation  of  the  H2O  fundamental  vibrational  modes  are  shown 
using  the  CBBB  surfaces.  Several  results  are  shown  including  those  shown  previously  using  the  BRSG 
surface  and  Gaussian  binned.  For  the  CBBB  results,  the  statistical  error  in  the  cross  sections  (1  a)  is 
approximately  13%  at  l.xlO'  m  and  41%  at  l.xlO'  m.  The  experimental  results  of  Dunn  et  a/,  are 
shown  as  well  as  the  QCT  results  of  Redmon  and  Schatz^  and  the  approximate  quantum  scattering 
calculations  of  Johnson,^  both  of  which  utilize  the  BRSG  surface. 

Examining  the  bending  excitation  cross  sections,  the  present  QCT  ZPO  results  are  un -physically 
large  and  remain  at  a  near-constant  value  over  the  entire  velocity  range.  As  with  the  BRSG  ZPO  result, 
these  large  cross  sections  are  due  to  borrowing  zero-point  energy  from  the  stretching  modes.  Maintaining 
zero-point  energy  (ZPl  and  ZP2  curves  which  we  have  suppressed  for  clarity)  considerably  lowers  the 
cross  section.  Likewise,  the  GB  approach  corrects  for  this  zero-point  energy  borrowing  deficiency.  At  6 
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and  8  km  s  ’,  the  CBBB  GB  and  BRSG  GB  bending  excitation  cross  sections  agree  well.  Closer  to 
threshold,  the  BRSG  GB  cross  section  is  a  factor  of  3-5  larger  than  the  CBBB  GB  cross  section,  but  the 
statistics  here  are  relatively  poor.  Both  the  CBBB  GB  and  BRSG  GB  cross  sections  agree  fairly  well  with 
the  approximate  quantum  scattering  results  of  Johnson,^  except  at  2  km  s  ’,  where  the  classical  results 
over-estimate  the  cross  section.  The  present  GB  results  with  both  surfaces  and  the  approximate  quantum 
scattering  results  are  all  about  10  times  smaller  than  the  measurement  of  Dunn  et  al.^ 

Figure  6  also  shows  results  for  the  stretching  vibrations.  Here  the  picture  is  more  complex.  The 
correction  in  going  from  the  ZPO  to  GB  result  for  the  symmetric  stretch  cross  section  is  large,  and  for  the 
asymmetric  stretch  it  is  very  small.  We  found  that  the  large  ZPO  symmetric  stretch  cross  section  is  due  to 
a  just  a  few  trajectories  on  the  S3  (highest  energy)  CBBB  surface.  The  S3  surface  is  the  only  one  that 
does  not  contain  an  H  -i-  OOH  channel,  and  symmetric  stretching  may  be  sensitive  to  this  difference. 
Gaussian  binning  completely  suppresses  these  trajectories.  However,  the  statistical  uncertainties  for  the 
symmetric  stretch  CBBB  result  are  large.  For  the  CBBB  GB  result,  the  asymmetric  stretch  cross  section  is 
much  larger  than  the  symmetric  stretch.  For  the  BRSG  surface  the  symmetric  and  asymmetric  stretch 
cross  sections  are  comparable.  For  both  the  stretches,  the  difference  between  the  CBBB  and  BRSG 
surface  cross  sections  are  about  a  factor  of  4  or  5.  The  stretch  excitations  are  more  sensitive  than  bending 
to  differences  in  the  non-reactive  versus  reactive  potential,  as  these  modes  sample  the  breaking  and 
formation  of  bonds  during  reaction.  The  measurement  of  Dunn  et  al.^  for  the  (001)  cross  section  at  4.3  km 
s  ’  is  at  least  several  orders  of  magnitude  above  the  present  results. 

Figure  7a  shows  H2O  vibrational  distributions.  The  peak  of  the  distributions  occurs  at  the  H2O 
zero  point  energy.  Vibrational  energy  in  excess  of  30  kcal  mol  ’  is  observed  at  high  collision  velocity  (8 
and  10  km  s’).  Figure  7b  shows  the  state -resolved  vibrational  cross  sections  as  a  function  of  the 
vibrational  energy.  At  the  highest  collision  velocities  an  exponential  fall-off  of  the  cross  section  with 
vibrational  energy  is  observed.  Figure  7c  shows  the  differential  angular  cross  sections  for  the  ground  and 
fundamental  vibrational  modes  at  8  km  s  ’  using  the  GB  method.  For  the  ground  state,  the  H2O  is  strongly 
forward  scattered  and  is  mostly  backscattered  for  the  fundamental  vibrations.  The  results  for  the  (100) 
mode  are  not  statistically  significant  and  so  we  do  not  show  them  here. 

In  Figure  8a,  the  rotational  cross  sections  for  H2O  are  shown  as  a  function  of  the  H2O  rotational 
angular  momentum,  J,  using  the  CBBB  potential.  These  have  been  summed  over  all  vibrational  states. 
QCT  ACT  GB  results  are  shown.  Similar  to  the  results  of  the  BRSG  potential,  as  the  collision  velocity 
increases,  the  higher  rotational  states  are  populated,  with  significant  rotational  cross  sections  above  J=20 
at  6,  8  and  10  km  s  ’.  At  collision  velocities  of  2  and  4  km  s  ’,  the  distribution  of  the  J  states  differs 
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slightly  from  the  BRSG  calculated  distrihutions.  With  the  CBBB  surfaces,  there  are  more  states  populated 
between  J=5  and  J=10.  Slight  differences  in  the  potential  surfaces  could  he  the  source  of  such  differences, 
especially  for  these  low  energy  excitations  where  the  rotational  energies  for  J<10  are  less  than  5  kcal 
mol'\  In  Figure  8h,  the  ro-vihrationally  resolved  cross  sections  are  given  for  the  ground  and  (010)  and 
(001)  vibrations  at  8  km  s  ’  using  the  QCT  ACT  GB  method.  The  (000)  cross  section  decreases 
monotonically  with  J,  but  the  (010)  cross  section  peaks  near  J=10,  while  the  (001)  cross  sections  are 
relatively  flat  with  J.  Figure  8c  shows  the  rotational  energy  distributions  for  the  ground  and  the  (010)  and 
(001)  modes.  Except  for  the  lowest  energy  bin,  these  show  an  exponential  fall-off  with  energy  consistent 
with  an  energy  gap  law. 

B.  Reactive  cross  sections 

In  Figure  9a,  the  reactive  cross  sections  are  shown  for  processes  in  Eqs.  1-3:  OH  formation,  H 
elimination,  and  O  atom  exchange.  We  have  discussed  these  cross  sections  in  detail  in  Briefly,  the 
OH  -I-  OH  reaction  dominates  the  reactive  cross  section  at  low  collision  velocities.  When  the  collision 
energy  increases  above  the  activation  barrier  for  H  elimination  at  8  km  s  ’,  the  cross  section  for  the  H  -i- 
OOH  reaction  becomes  significant  and  competitive  with  the  hydrogen  abstraction  reaction  at  10  km  s  ’. 
The  0-atom  exchange  process  has  a  relatively  low  barrier  (33  kcal  mol  ’)  and  turns  on  at  6  km  s  ’.  The 
opacity  function  is  given  for  each  of  the  reactions  at  a  collision  velocity  of  10  km  s  ’  in  Eigure  9b.  The 

opacity  function,  P(b),  has  the  following  relation  to  the  cross  section,  cr,  :  cr  =  P(h^hdh,  where  b  is 
the  impact  parameter  and  b^ax  is  the  maximum  impact  parameter.  Low  impact  parameters  correspond  to 
0-atom  interaction  in  close  proximity  to  the  H2O  center  of  mass.  With  sufficient  energy  at  8  km  s  ’,  these 
low  impact  parameter  interactions  tend  to  lead  to  the  H  -1-  OOH  product  channel.  As  the  impact  parameter 
decreases,  the  probability  of  the  H  -1-  OOH  channel  increases  indicating  a  ‘direct’  interaction  is  necessary. 
The  hydrogen  abstraction  reaction,  however,  occurs  at  larger  impact  parameters;  the  formation  of  OH  -1- 
OH  dominates  the  cross  section  at  impact  parameters  greater  than  ~1  au.  The  opacity  function  for  the  O- 
atom  exchange  reaction  is  similar  to  that  of  the  H  -1-  OOH  reaction  albeit  two  orders  of  magnitude  lower. 
This  indicates  that  exchange  reaction  requires  a  low  impact  parameter  between  the  O  and  H2O  and  in 
addition  is  limited  to  particular  geometries.^  Eor  these  total  reaction  cross  sections,  differences  between 
ZPO,  ZPl,  and  ZP2,  and  GB  results  were  very  small. 

We  now  examine  the  ro-vibrationally  resolved  cross  sections  and  product  energy  distributions  of 
the  OH  formation  and  H  elimination  reactive  channels.  Statistics  for  the  ro-vibrationally  resolved  0-atom 
exchange  channel  were  too  poor  to  present  results  for  this  channel.  Results  for  QCT  ACT  with  Gaussian 
binning  (GB)  will  be  shown.  Very  near  threshold  the  ZPO  results  produced  small  cross  section  products 
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that  were  only  energetically  allowed  by  borrowing  zero-point  energy.  Besides  this,  differences  between 
ZPO,  ZPl,  ZP2,  and  GB  results  were  generally  small. 

1.  OH-i-OH  reaction 

The  vibrational  energy  distributions  for  the  OH  -i-  OH  products  are  shown  in  Figure  10a.  The 
contributions  for  each  OH  product  has  been  summed  and  divided  by  two,  so  the  integrated  distribution  at 
a  given  collision  velocity  equals  the  reactive  cross  section.  The  distributions  peak  near  the  zero  point  of 
OH  near  ~5  kcal  mol'^  At  the  highest  collision  velocities  examined,  cross  sections  are  observed  at 
vibrational  energies  in  excess  of  20  kcal  mol  '.  Beyond  the  zero-point  peak,  this  energy  distribution  falls 
off  approximately  exponentially.  This  is  more  clearly  seen  in  Figure  10b,  which  shows  the  vibrationally 
resolved  cross  section  versus  vibrational  energy,  where  again  the  two  product  OH  cross  sections  are 
summed  and  divided  by  two.  In  a  previous  study, we  have  shown  that  two  distinct  types  of  OH  products 
are  formed.  A  spectator  OH,  corresponding  to  the  OH  H2O  bond  that  is  not  broken,  occurs  with  lower 
vibrational  and  rotational  energy,  and  an  active  OH,  corresponding  to  the  OH  bond  that  is  newly  created, 
is  formed  with  higher  vibrational  and  rotational  energy.  This  bi-modal  behavior  can  be  seen  in  the  Figure 
10b,  where  after  the  OH(v=0)  peak,  there  is  an  exponential  fall-off  of  the  vibrational  energy.  However, 
the  OH(v=0)  level  rises  above  this  exponential  line  for  a  given  collision  velocity.  This  is  due  to  the 
contribution  of  the  OH  spectator  which  is  formed  mostly  in  the  v=0  state,  while  the  active  OH  is  formed 
with  a  wide  distribution  of  vibrational  levels. 

The  rotationally  resolved  cross  sections  for  OH  are  shown  in  Figure  1 1  where  contributions  from 
both  OH  products  have  been  summed  and  divided  by  2.  In  Figure  11a  we  show  the  J-resolved 
distributions  summed  over  vibrational  levels.  At  the  higher  collision  velocities,  these  distributions  are 
actually  sums  of  a  spectator  OH  distribution  peaked  at  low  J  (between  J=5  and  10)  and  an  active  OH 
distribution  which  is  broad  with  a  higher  J  tail.^^  This  can  be  seen  in  a  different  way  in  Figure  11b,  which 
shows  the  J  resolved  cross  sections  at  8  km  s  ’  for  the  OH(v=0)  and  OH(v=l)  levels.  For  OH(v=0),  the 
peak  near  J=8  is  due  to  the  spectator  OH  and  the  higher  energy  J  tail  is  due  to  the  active  OH.  For 
OH(v=l),  the  J-resolved  cross  section  has  poor  statistics.  But  the  very  broad  distribution  versus  J  is 
evident  and  is  due  mostly  to  active  OH  formed  in  the  v=l  level.  Figure  11c  shows  the  rotational  energy 
distributions  at  8  km  s  ',  for  OH(v=0, 1,2,3).  At  the  highest  collision  velocity,  these  distributions  extend 
beyond  20  kcal  mol  '.  Above  5  kcal  mol  '  they  can  be  approximately  described  by  an  exponential  fall  off 
with  energy.  This  view  of  the  cross  section  in  Figure  11c  is  basically  a  remapping  for  Figure  11b,  where 
the  low  J  values  have  been  compressed  to  the  lower  energy  bin  and  the  higher  J  values  are  spread  across 
widely  to  higher  energy  bins. 
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2.  OOH  +  H  reaction 


Figure  12  shows  the  vihrationally  resolved  cross  sections  for  the  OOH  product.  To  our 
knowledge  this  is  the  first  time  vihrationally  resolved  cross  sections  have  been  reported  for  this  channel. 
Figure  12  shows  the  vibrational  energy  distributions  for  OOH  at  a  collision  velocity  of  10  km  s  ’,  where 
we  show  the  separate  contributions  from  the  X  state  of  OOH  (energy  threshold  of  ~  60  kcal  mol  ')^  and 
the  A  state  of  OOH  (energy  threshold  of  ~84  kcal  mol  ').^  The  zero-point  energy  of  these  states  is  about 
~10  kcal  mol  '.  In  contrast  to  the  vibrational  distributions  of  H2O  and  OH  which  are  strongly  peaked  near 
their  zero-point,  the  vibrational  energy  distribution  for  OOH  is  extremely  broad.  The  X  state,  having  a 
lower  energy  threshold,  dominates  the  distribution  and  extends  up  to  40  kcal  mol  ’.  We  have  noted 
previously  that  the  H  elimination  channel  is  unusual  in  that  -40%  of  available  collision  energy  ends  up  in 
the  internal  energy  of  the  OOH  product.^  The  vihrationally  resolved  cross  section  is  shown  in  Figure  12b, 
where  the  X  state  (which  dominates)  and  A  state  contributions  have  been  summed.  The  harmonic  normal 
modes  for  the  X  state  of  OOH  have  energies  of  3863.7(H-stretch),  1666.7(bend),  and  1041.8  (0-0 
stretch)  cm  ’  and  we  use  the  notation  from  high  to  low  energy  (H-stretch,  bend,  0-0  stretch).  The  cross 
sections  for  the  ground  (000)  and  first  excited  bend  (010)  (which  is  not  the  lowest  energy  mode)  modes 
are  nearly  equal,  which  is  consistent  with  the  broad  vibrational  energy  distribution.  Below  values  of  1x10 
m^  (which  are  not  shown),  there  is  excitation  of  a  large  number  of  states  at  higher  vibrational  energies, 
which  make  up  the  high  energy  tail  of  the  X  state  vibrational  energy  distributions.  These  states  have  large 
statistical  uncertainties  and  require  many  more  trajectories  to  give  a  quantitative  estimate  of  their  cross 
sections.  The  angular  distribution  cross  sections  for  the  ground  (000)  and  (001)  and  (010)  states  at  10  km 
s  ’  are  shown  in  Figure  12c.  The  angle  is  the  center  of  mass  scattering  angle.  High  scattering  angles 
correspond  to  scattering  of  the  OOH  product  in  the  same  direction  as  the  reactant  H2O  molecule.  The 
dominant  backscattering  component  for  this  channel  has  been  seen  previously.’’^  Reaction  takes  place  at 
small  impact  parameters  (hard  collisions). 

In  Figure  13a,  we  show  the  J-resolved  cross  section  for  the  OOH  product  at  10  km  s  ’.  The  cross 
section  is  summed  over  X  and  A  states  and  all  vibrational  states.  The  J-resolved  cross  section  has  a  broad 
distribution  with  significant  populations  from  J=20  to  J=70  and  extends  past  J=100.  The  broad 
distribution  extending  to  high  J  is  due  to  two  factors,  the  relatively  large  amount  of  energy  deposited  in 
internal  motion  and  the  rotational  constants  of  OOH.  The  experimentally  determined  rotational  constants 
for  OOH  (X  state)  are  A=20.35,  B=1.12,  and  C=1.06  cm  ’.^®  This  corresponds  to  a  near  prolate 
symmetric  top,  with  A  >  B  =  C  and  energy  levels  E^o^j  =  BJ(J  -I-  1)  -I-  (A  —  B)K^,  -J. 

Because  of  the  relatively  small  value  of  the  smallest  rotational  constants,  for  K=0  very  large  J  values  can 
be  populated.  Also,  because  of  the  large  difference  between  the  largest  and  smallest  rotational  constants,  a 
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given  J  level  can  give  rise  to  a  large  number  of  K  states  which  span  a  large  energy  range.  This  is  shown  in 
Figure  13c,  where  rotational  energy  is  plotted  as  a  function  of  J.  For  rotational  energies  of  30  kcal  mol  ' 
and  lower,  possible  angular  momentum  states  are  possible  up  to  J=100  for  small  K  values.  Figure  13b 
shows  the  rotational  energy  distribution  summed  over  the  X  and  A  states  and  all  vibrational  states,  which 
is  nearly  flat  with  rotational  energy.  This  distribution  is  fundamentally  different  than  the  corresponding 
H2O  or  OH  distributions  where  the  low  energy  population  (<  10  kcal  mol  ')  is  at  least  an  order  of 
magnitude  greater  than  the  high  energy  population  (>  10  kcal  mol  '). 

3.  Rate  constant  calculations 

As  a  test  on  the  CBBB  potential.  Figure  14  shows  the  computed  rate  constants  versus  temperature 
using  the  CBBB  surfaces  for  the  reactions  of  (a)  O  +  H2O  OH  +  OH  and  (b)  OH  +  OH  H2O  +  O 
along  with  measured  data  For  the  O  +  H2O  rate  constant,  results  using  an  older  set  of  potential 
surfaces^^  using  QCT  are  also  shown.  For  the  present  CBBB  results  for  the  O  +  H2O  rate  constant, 
contributions  from  each  of  the  three  electronic  states  are  summed  and  divided  by  3.  The  CBBB  result 
agrees  well  with  the  highest  valued  experimental  measurement.  For  the  OH  +  OH  rate  constant, 
contributions  from  each  electronic  state  are  multiplied  by  1/4  to  reflect  the  fact  the  OH  +  OH  reactants  are 
4x  degenerate.  Our  calculations  neglect  contributions  from  the  fourth  electronic  state,  but  we  have  seen 
that  the  contribution  from  the  3rd  highest  state  is  negligible,  so  we  expect  that  the  fourth  state  we  are 
leaving  out  does  not  contribute  significantly.  The  present  results  agree  fairly  well  with  the  highest  valued 
measurements. 

V.  CONCLUSIONS 

Classical  dynamics  calculations  for  0(^P)  +  H20('Ai)  hyperthermal  collisions  were  performed. 
Several  methods  were  used  to  produce  ro-vibrationally  state -resolved  product  cross  sections  and  to 
enforce  zero-point  maintenance  from  analysis  of  the  classical  trajectories.  In  one  series  of  calculations,  a 
non-reactive  surface  (BRSG)  ’  ’  used  in  past  classical  and  quantum  collision  studies  was  used.  Results 
with  this  surface  for  collisional  excitation  suggest  that  classical  dynamics  calculations  which  include 
Gaussian  binning  procedures  and/or  selected  zero-point  maintenance  algorithms  can  produce  results 
which  approximate  quantum  scattering  cross  sections.  Without  these  procedures,  the  classical  cross 
sections  can  be  many  orders  of  magnitude  greater  than  the  quantum  cross  sections  for  exciting  the 
bending  vibration  of  H2O,  especially  near  threshold.  This  is  due  to  energy  borrowing  from  stretching 
modes  which  dip  below  zero-point  values.  In  a  second  series  of  calculations,  recently  developed  reactive 
0(^P)  -I-  H20('Ai)  surfaces  were  used  (CBBB).^  For  these  surfaces  the  three  lowest  triplet  states  of  the  O 
-I-  H2O  system  included  reaction  to  form  OH  -1-  OH,  H  -1-  OOH,  and  0-atom  exchange.  For  bending 
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vibrational  excitation,  results  for  the  BRSG  and  CBBB  surfaces  agreed  fairly  well.  For  the  stretches  the 
cross  sections  showed  differences  of  factors  of  4  or  5,  although  the  statistical  uncertainties  were  large. 
The  stretch  excitations  are  more  sensitive  than  bending  to  differences  in  the  non-reactive  versus  reactive 
surfaces,  as  these  modes  sample  the  breaking  and  formation  of  bonds  during  reaction.  The  measurements 
of  Dunn  et  al.^  for  the  (010)  and  (001)  cross  sections  at  4.3  km  s  ’  were  at  least  an  order  of  magnitude 
above  the  present  results.  For  the  reactive  channels,  the  present  calculations  using  the  CBBB  surfaces 
show  that  at  higher  velocities  there  is  an  unusually  large  amount  of  product  internal  excitation.  For  OOFl, 
where  40%  of  available  collision  energy  goes  into  internal  motion,  the  excited  product  vibrational  and 
rotational  energy  distributions  are  relatively  flat.  Other  product  channel  distributions  show  an  exponential 
fall-off  with  energy  consistent  with  an  energy  gap  law.  The  present  detailed  distributions  and  cross 
sections  can  serve  as  a  guide  for  future  hyperthermal  measurements  of  this  system. 
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TABLES 


Table  I.  Parameters  used  in  QCT  calculations. 


Parameter 

Value 

Collision  velocities 

2,  4,  6,  8  (BRSG)  2,  4,  6,8,  10  (CBBB)  km  s  ' 

Maximum  impact  parameter 

7.54  au,  atot=  1.78.60  au  =  5.03x10'’*’  m^ 

Initial  separation 

17.0  au 

Number  of  trajectories  per  velocity 

UxlO^CBRSC),  l.xlO^  (CBBB  per  surface) 

Statistical  Uncertainty  (la) 

10%  at  5.xl0 ""  m’’;  32%  at  5.xl0 m’’  (BRSG) 

13%  at  l.xlO'^^  m^;  41%  at  l.xlO'^^  m^  (CBBB) 

Time-step 

1x10  s 

Product  analysis  time 

5.12x10  ’'' s 

Maximum  trajectory  time 

l.xlO”  s 

Rotational  energy  bin  width 

2  kcal  mof’ 

Translational  energy  bin  width 

2  kcal  mof’ 

Vibrational  energy  bin  width 

2  kcal  mof’ 

Angular  bin  width 

10  degrees 

z.p.  lower  limit 

50  cm  ’ 

z.p.  dV/dx^  (energy/length^)  threshold 

3.16x10  '  au  (BRSG),  3.16x10”  au  (CBBB) 

Average  energy  conservation  (8  km  s  ’) 

0.19  kcal  mof’  (BRSG)  6.8x10”  kcal  mof ’(CBBB) 
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Table  II.  Frequencies  (cm  *)  and  equilibrium  geometries  (au,  radians)  of  reactants  and  products.  The 
BRSG  and  CBBB  frequencies  are  harmonic  frequencies  determined  from  these  surfaces.  The  literature 
value  frequencies  are  derived  from  measurements  for  excitation  of  fundamental  modes  and  therefore 
include  anharmonic  and  higher  order  contributions. 


Molecule 

BRSG  and  CBBB  Surfaces 

Literature 

H2O 

Freq.=3863.7,  1688.1,3973.2  (BRSG) 

Roh=  1.81,  0HOH  =1.83  (BRSG) 

Freq.=  3753.5  ,1781.4,3760.5  (CBBB)  SI 

Roh=1.83,  0hoh=1.81  (CBBB)  SI 

Freq.=3657,  1595,  3756 "" 

Roh=1.81 , 0HOH  =1.82 

OOH 

Freq.=  3863.7,  1666.7,  1041.8  (CBBB)  SI  (X) 

Roo=2.56,  Roh=1.83,  Oqoh  =1.77 

Freq.=3436,  1392,  1098 

Roo  =2.51,  Roh  =1.8,  OooH  =1.82'' 

OH 

Freq.  =  3859.9  (CBBB)  SI 

Roh=  1.85 

Freq.=3738 

Ron  =1.83"“ 
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FIGURE  CAPTIONS 


Figure  1.  Potential  energies  for  0(^P)  +  H20(^Ai)  using  the  BSRG  surface  of  Refs.  7,  8,  and  18  and  the 
CBBB  SI  surface  of  Ref.  9.  In  (a)  and  (h)  H2O  is  held  near  its  equilihrium  geometry  Roh=1-83  au,  0hoh 
=1.82  rad  and  the  0-atom  is  moved  in  the  H2O  plane.  In  (c)  and  (d)  the  H2O  geometry  corresponds  to 
bond  lengths  of  RoHi=l-83au  and  Roh2=2.40  au,  and  an  angle  0hoh  =1.82  rad.  The  0-atom  is  moved  in 
the  H2O  plane. 

Figure  2.  Fundamental  vibrational  excitation  cross  sections  for  O  -1-  H2O  collisions  using  the  BRSG 
surface.  Please  see  text  for  description  of  data. 

Figure:  3  Contour  plots  of  the  logl0  of  the  vibrational  excitation  probabilities  normalized  to  1.0  for  each 
sub-panel  versus  the  bending  and  averaged  stretching  vibrational  actions  using  the  BSRG  surface. 

Figure  4:  Vibrational  cross  sections  for  O  -1-  H2O  collisions  with  the  BSRG  surface:  (a)  Vibrational 
energy  distribution  cross  sections  (m^/kcal  mof^)  at  2,  4,  6,  and  8  km  s  *  for  QCT  ZPl.  (b)  Vibrational 
excitation  cross  sections  at  8  km  s"^  and  6  km  s'^  for  QCT  GB.  (c)  Angular  distribution  cross  sections  for 
the  ground  and  fundamental  H2O  modes  at  8  km  s‘^  for  QCT  GB. 

Figure  5:  Rotational  cross  sections  for  O  -1-  H2O  collisions  with  the  BSRG  surface,  (a)  J-resolved  cross 
sections  versus  collision  energy;  below  4  kcal  mol-1,  results  are  shown  from  Ref.  28,  (b)  J-resolved  cross 
sections  at  a  number  of  collision  velocities,  (c)  J-resolved  cross  sections  for  the  ground  and  vibrational 
fundamental  modes  at  8  km  s  ',  (d)  Rotational  energy  distributions  for  the  ground  and  vibrational 
fundamental  modes  at  8  km  s  ’. 

Figure  6:  Vibrational  excitation  cross  sections  for  O  -1-  H2O  collisions,  (a)  excitation  to  the  (010)  state  (b) 
excitation  to  the  (100)  state,  (c)  excitation  to  the  (001)  state.  Black  line  with  circles  (CBBB  GB)  denotes 
Gaussian  binned  results  using  the  CBBB  surfaces;  Red  line  with  squares  (CBBB  ZP0)  denotes 
conventional  histogram  and  no  zero-point  algorithm  using  the  CBBB  surfaces;  Blue  line  with  triangles 
(BRSG  GB)  denote  present  QCT  Gaussian  binned  results  using  the  BRSG  surface;  Purple  lines  denote 
results  from  Ref.  7;  Purple  lines  with  symbols  are  tabulated;  Purple  lines  without  symbols  are  a  function 
fit;  Green  lines  with  diamonds  (VCC-IOS)  are  the  approximate  scattering  results  Johnson®;  Cyan  blue 
symbols  (Dunn  et  al.)  are  shock  tube  measurements  from  Dunn  et  al.^ 

Figure  7:  Fundamental  vibrational  cross  sections  for  O  -1-  H2O  collisions  with  the  CBBB  surfaces:  (a) 
Vibrational  energy  distribution  cross  sections  (m^/kcal  mol  ')  at  2,  4,  6,  8,  and  10  km  s  '  for  QCT  ZPl.  (b) 


29 


Distribution  A:  Approved  for  pubiic  reiease;  distribution  uniimited. 


Vibrational  excitation  cross  sections  at  6,  8,  and  10  km  s  ’  for  QCT  ACT  GB.  (c)  Angular  distribution 
cross  sections  for  the  (000),  (010),  and  (001)  H2O  modes  at  8  km  s  ’  for  QCT  ACT  GB. 

Figure  8:  Rotational  cross  sections  for  O  +  H2O  collisions  with  the  CBBB  surfaces  leading  to  excited 
H2O.  (1)  J-resolved  cross  sections  at  a  number  of  collision  velocities,  (b)  J-resolved  cross  sections  for  the 
(000),  (010),  and  (001)  modes  at  8  km  s  ’,  (c)  Rotational  energy  distributions  for  the  (000),  (010),  and 
(001)  modes  at  8  km  s  ’. 

Figure  9:  (a)  The  total  reactive  cross  sections  versus  collision  velocity  (b)  The  opacity  function  versus 
impact  parameter  at  10  km  s  ’  for  each  reactive  channel. 

Figure  10:  Vibrational  cross  sections  for  formation  of  OH  +  OH  using  the  CBBB  surfaces,  (a) 
Vibrational  energy  distributions,  (b)  Vibrationally  resolved  cross  sections. 

Figure  11:  Rotational  cross  sections  for  O  +  H2O  collisions  with  the  CBBB  surfaces  leading  to  OH  +  OH. 
(a)  J  resolved  cross  sections  summed  over  vibrational  states,  (b)  The  ro-vibrationally  resolved  cross 
section  ground  and  first  excited  states  at  8  km  s  ’,  (c)  Rotational  energy  distributions  for  OH(v=0,l, 2,3,4) 
at  8  km  s  ’. 

Figure  12:  Vibrational  cross  sections  for  O  +  H2O  collisions  with  the  CBBB  surfaces  leading  to  H  + 
OOH.  (a)  Vibrational  energy  distributions  for  the  X  and  A  states  of  OOH  at  10  km  s  ’,  (b)  Vibrationally 
resolved  cross  sections  10  km  s  ’,  (c)  The  differential  angular  cross  section  for  the  (000),  (001),  and  (010) 
states  at  10  km  s  ’. 

Figure  13:  Rotational  cross  sections  for  O  +  H2O  collisions  with  the  CBBB  surfaces  leading  to  H  +  OOH. 
(a)  J-resolved  cross  section  for  OOH  at  10  km  s  ’  (b)  Rotational  energy  distributions  of  OOH  at  10  km  s  ’, 
(c)  Rotational  energy  of  OOH  as  a  prolate  symmetric  top  versus  J. 

Figure  14:  (a)  Rate  constants  for  the  reaction  of  O  -1-  H2O  OH  -1-  OH  in  m^  molecule  ’  s  ’  as  a  function 
of  temperature  in  Kelvin.  The  legend  notation  is  as  follow:  JCP2010  are  the  present  CBBB  results. 
JCP2005  are  the  results  of  Ref.  35.  LIF  1991  is  from  Ref.  21.  TSA  1986  is  from  Ref.  22.  (b)  Rate 
constants  for  the  reaction  of  OH  -1-  OH  ^  O  -1-  H2O  in  m^  molecule  ’  s  ’  as  a  function  of  temperature  in 
Kelvin.  JCP  2010  are  the  present  CBBB  results.  WOO  1995  is  from  Ref.  23.  BAU  1992  is  from  Ref.  20. 
TSA  1986  is  from  Ref.  22. 
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Figure  2 
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Figure  3 
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Figure  4 
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Figure  6 
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Figure  7 
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Figure  8 
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Cross  Section  (m^) 


Figure  9 
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Figure  10 
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Figure  12 
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Figure  13 
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Figure  14 
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